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We report on the manipulation of the center-of-mass motion ('sloshing') of a Bose Einstein con- 
densate in a time-averaged orbiting potential (TOP) trap. We start with a condensate at rest in 
the center of a static trapping potential. When suddenly replacing the static trap with a TOP trap 
centered about the same position, the condensate starts to slosh with an amplitude much larger 
than the TOP micromotion. We show, both theoretically and experimentally, that the direction of 
sloshing is related to the initial phase of the rotating magnetic field of the TOP. We show further 
that the sloshing can be quenched by applying a carefully timed and sized jump in the phase of the 
rotating field. 
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I. INTRODUCTION 

Time averaged potentials (TAP) offer a versatile tool 
for trapping both charged and neutral particles. For 
neutral atoms the most common example in this class 
of traps is the Time-averaged Orbiting Potential (TOP) 
which was used in the experiments in which the first 
Bose-Einstein condensate (BEC) was created [H, 0). The 
TOP trap consists a magnetic quadrupole trap [3, H| 
shifted by a uniform magnetic modulation field rotating 
at a high (audio) frequency. As this rotation is slow as 
compared to the Larmor precession of the atomic mag- 
netic moments, the atoms remain polarized with respect 
to the instantaneous effective magnetic field [j| as follows 
from the adiabatic theorem. On the other hand the ro- 
tation is fast as compared to the orbital motion of the 
atoms. As a consequence, the atomic motion consists of 
a fast rotating part (micromotion), superimposed on a 
slow oscillating part (macromotion). In the simplest the- 
oretical description, the static approximation, the micro- 
motion is eliminated by time-averaging the instantaneous 
potential over a full cycle of the modulation field. 

Suppose we load a particle with given momentum p at 
position ro in a TOP trap using a sudden switch-on pro- 
cedure. One might naively guess that the ensuing motion 
is given by the dynamics in the time-averaged potential, 
subject to the initial conditions r = ro and p = Po but 
this guess turns out to be wrong. In fact, one can show 
that the initial conditions for the slow motion depend on 
the phase of the TOP at the time of switch-on. This phe- 
nomenon was analyzed by Ridinger and coworkers [y, 0] 
for the special case of a one-dimensional rapidly oscil- 
lating potential (ROP) with zero average. Ridinger et 
al. also showed, first for a classical particle @ and subse- 
quently for the quantum case Q , that the amplitude and 
energy associated with the slow motion can be altered by 
applying a suitable phase jump in the rapidly oscillating 
field. 

In this paper we show, both theoretically and exper- 
imentally, that the dependence on initial phase and the 
possibility to influence the motion by phase jumps, is also 



present for a two-dimensional rotating TOP field. In par- 
ticular we show that a cloud of atoms which is initially 
at rest with zero momentum acquires a sloshing motion 
as soon as the TOP is suddenly switched on. This is true 
even if the cloud is initially at the minimum of the effec- 
tive potential. The amplitude of this slow macromotion 
is much larger than that of the fast micromotion while 
the direction of sloshing depends on the TOP phase at 
switch-on. We also demonstrate that this macromotion 
can be almost entirely quenched by applying a carefully 
timed and sized phase jump in the TOP field. 

The motion of atoms and ultracold atomic clouds in 
TOP traps have been extensively described in the liter- 
ature. Following the achievement of the first BEC Q, 
the use of the axially symmetric TOP was described 



theoretically in 
other groups 



U3| and explored experimentally by 
to study properties of the BEC. The 
idea of the TOP was extended to an asymmetric triax- 
ial TOP trap developed by and also used by other 
groups [13, EH- Further a number of other variations 
were introduced: In many cases, it turns out to be con- 
venient to switch on the TOP after a preparative stage 
of cooling in a conventional static trap such as a mag- 
netic quadrupole trap (see e.g. [HI), an optically plu gged 
magnetic quadrupole [2(| and Ioffe-configurations [2l|— 
[23j . Often, the transfer of the cloud from the static to the 
TOP trap cannot be performed adiabatically for topolog- 
ical reasons. Bearing this in mind, it becomes relevant to 
carefully analyze the dynamics that may be induced by 
a sudden switch-on of the TOP. In addition, applications 
which require manipulation of a BEC are heavily depen- 
dent on precise control of the location of the atomic cloud 
and can thus benefit from the techniques described. 

In our experiments the condensate is prepared in a 
Ioffe-Pritchard (IP) trap before transferring to a TOP 
because the use of radio-frequency (rf ) induced evapora- 
tive cooling is more efficient in a static magnetic trap, 
resulting in larger condensates. Once transferred to the 
TOP we can create trapping geometries that are difficult 
to realize using a static magnetic potential without intro- 
ducing Majorana losses associated with the presence of 
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zero-field points. An example is the double well potential 
used in [2l| . 

The remainder of this paper is organized as follows. In 
Section |TT] we calculate the motion of a cloud of atoms 
in a TOP which at switch-on is at rest at the center of 
the trap. We discuss the motion that results and de- 
rive the conditions under which a phase jump can lead 
to a substantial reduction of the energy associated with 
the slow motion of the cloud. In Section Mil we discuss 
the experimental details and the preparation of the BEC 
and its transfer to the TOP. In Section IIVI we present 
the experimental results and compare with the theory of 
Section [TXJ Finally in Section [V] we give a summary and 
conclusion. 



II. THEORY 
A. Time-averaged Ioffe-Pritchard potential 

In the literature the term TOP is most often used for a 
spherical-quadrupole trap combined with a rotating uni- 
form magnetic modulation field. In this paper we will 
use the term TOP in a broader context, to include the 
magnetic trapping potential created by combining a IP 
trap with rotating modulation field. Challis et al. [l3j 
have shown that the dynamical eigenstates of a degener- 
ate Bose gas in a TOP are given by solutions of the usual 
Gross-Pitaevskii equation but taken in a circularly trans- 
lating reference frame, that is, a reference frame the ori- 
gin of which performs a rapid circular motion but retains 
a constant orientation. In particular this implies that the 
center of mass of a condensate in its ground state per- 
forms the same micromotion in a TOP as a point particle 
with the magnetic moment of an atom. In this spirit we 
use as a 87 Rb condensate to study the micromotion and 
macromotion in a TOP. 

We consider a cigar-shaped Ioffe-Pritchard potential 

U(q, z) = /^aV + (Bo + !/3* 2 ) 2 , (1) 

where g(t) is the radial position of a test atom with re- 
spect to the IP symmetry axis, p, the magnetic moment 
of the atom, and a, f3, Bq the parameters for the radial 
gradient, the axial curvature and offset value of the IP 
magnetic field. Eq. (TTJ) represents an approximate expres- 
sion for the IP trap which is valid for a 2 ^> (3Bq and in 
the limit g < a//3 |, [U H|. 

In the presence of the TOP field we transform to the 
circularly translating frame (l3j and have 

g(t) = {x - p m cos(W + (f) m ), y - p m sin(W + (f) m )}, (2) 

where {x, y, z} = {p, z} = r is the position of the atom 
in the laboratory frame and the IP symmetry axis is dis- 
placed over a distance p m — B m /a in the direction 

p m = {cos(ujt + (f) m ),sm(ujt + (j) m )} (3) 



by the uniform modulation field 

B m = B m {cos(iut + (j> m ), - sin(wt + cj> m )} (4) 

applied perpendicular to the z axis. The y axis is taken 
along the vertical direction, the xz plane being horizon- 
tal. The modulation field B m rotates at angular fre- 
quency — ui (phase —<f> m ) about the horizontal z axis as 
illustrated in Fig. [I] Notice that the sense of rotation 
of the IP-ficld-minimum is opposite to that of the B m 
field, in contrast to the original TOP configuration 
where the field-zero rotates in the same direction as the 
bias field. This reflects the difference between the 2D- 
quadrupole symmetry of the IP trap and the axial sym- 
metry of the spherical-quadrupole trap. The rotation of 
the modulation field B m also gives rise to a fictitious field 
B w which has to be added or subtracted from the offset 
field Bo, depending on the sense of rotation, 

B ^B (1±B W /B ) =B (1 ±u/u L ), (5) 

where ujl = gFpBBo/h is the Larmor frequency of mag- 
netic moment of the atoms, with gp the hyperfine g fac- 
tor and hb the Bohr magneton. In a standard TOP, 
the fictitious field in combination with gradient of the 
quadrupole field gives rise to a shift of the equilibrium 
position of the cloud in the direction of the axis around 
which the field rotates @, In our IP-TOP the axial 
field is homogeneous near the origin and the shift is ab- 
sent; the change in Bq turns out to be small and will be 
neglected in this paper. 

For (3 — and B = the potential U(g,z) corre- 
sponds to that of a two-dimensional quadrupole field with 
a zero-field line that rotates at distance p rn about the z 
axis as a result of the modulation. For Bo — the dis- 
tance p m is known as the radius of the 'circle of death'. 
For B < the potential corresponds to two TOP traps 
separated by Az = 2(2\B \/ P) 1 ^ 2 [2l|. In this paper we 
will consider only the case Bo > 0. 

In the common description of the TOP one analyzes 
the motion in an effective potential, obtained by time 
averaging the static trap over a full rotation period of the 
B m field. For Eq. (fT]) this procedure yields the effective 
potential 

i r** 

W(r) = —J U(x - p m cos C, y - p m sin £ z)d(, (6) 

where £ = ujt + (f> m . For the cigar-shaped IP potential we 
consider the condition 

lu > n p > n z , (7) 

where, for an atom of mass m, the quantity f2 z = 
(fx (3/m) 1 / 2 is the axial harmonic oscillation frequency in 
the effective potential U(0, 0, z). Analogously, harmonic 
oscillation frequency in the radial plane is given by 



fi, = W£ir(i-W/Bg) = n, ( 8 ) 
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where B = (B% + B^) 1 ' 2 is offset value of the effective 
potential at the origin [2l| . 

The first inequality in Eq. <[7|) ensures that the fast and 
slow radial motions of the atoms can be separated, which 
is the well-known operating regime for a TOP trap 
The second inequality implies that the axial motion in 
the effective trap is slowest and that the motion can be 
treated as quasi two-dimensional in the radial plane. 

To account for the acceleration due to gravity (g), the 
gravitational potential mgy has to be added to Eqs. (fTJ) 
and ([5]) ■ The main effect is to shift the minimum of the 
potentials in the negative y direction by the amount 



Ay = g/n 2 . 



(9) 



This expression holds as long as the gravitational sag Ay 
is much smaller than the harmonic radius pu = Bo /a. 

Since p% > p m , the effective potential © may be 
treated as harmonic as long as the motion is confined 
to a region around the z axis that is small compared 
to p rn . For our experiment the harmonic approximation 
holds rather well and is sufficient for gaining qualitative 
insight in the micro- and macromotion as will be shown 
in Section III Bl Refinements associated with switch-on 
transients and gravity are discussed in Appendix In 
the numerical analysis of Section III CI we solve the clas- 
sical equations of motion in the full time-dependent po- 
tential Eq. (TT]). In this context we also comment on the 
validity of the harmonic approximation. 
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FIG. 1: (color online) Schematic diagram of the magnetic field 
configuration in relation to the orbit of stationary micromo- 
tion (solid blue circle). The view is along the (horizontal) z 
axis. The orbital position and velocity of the micromotion are 
denoted by p — —poP m and Vq. The IP symmetry axis ro- 
tates at frequency uj (with initial phase <f> m ) about the z axis 
on the circle of radius p m (dashed black circle). Note that 
the TOP field B m = B m p m rotates at frequency — uj (phase 
—<f>m), reflecting the 2D-quadrupole symmetry (dashed red 
circle) of the IP trap. 



B. Micromotion and macromotion 

To analyze the effect of switching on the B m field at 
t = we first consider an atom 'at rest' in the center of 
the effective trapping potential U{p,z). Such an atom 
exhibits no period-averaged dynamics (no macromotion) 
but only circular micromotion at a frequency uj about the 
origin as illustrated in Fig.Q] The radius of this station- 
ary micromotion, 

/5o = S( 1+s ° 2/ ^" /2 ' (10) 

follows from the condition F c = muj 2 po for the centripetal 
force F c = -V P C/| P=0 = M«(l + B 2 / B^y 1 ' 2 p m . The 
speed of this stationary micromotion, 

Vo= ojp = ^(l + B 2 /B 2 m y 1/ \ (II) 
muj 

is directed orthogonally to the direction p m . Such pure 
micromotion only results if at t = the atom is already 
moving at speed vq along a circle of radius po about 
the origin and is located at position p = —poP m (see 
Fig.Q]). Obviously an atom at t = at rest at the origin, 
p = {0,0} does not satisfy these initial conditions and 
as a consequence its macromotion will start with a finite 
launch speed. We will see that the result is elliptical mo- 
tion at frequency f2, with the long axis approximately 



perpendicular to the initial direction of p m and with a 
substantial amplitude, of order (w/Q)p . Usually this 
motion is undesired and our aim is to quantify it and 
subsequently quench it by imparting a phase jump to 
the TOP-field. 

To gain insight into the way in which the sudden 
switch-on of the TOP influences the macromotion of an 
atom initially at rest at the origin, we first consider a 
simple model in which it is assumed that the motion in 
the radial plane can be decomposed into two harmonic 
components, oscillating at the micromotion and macro- 
motion frequencies u> and Q, respectively. The position 
pit) and velocity pit) are given by 

pit) = {po cos(wi + (f>), po sin(cj£ + 4>)} + 

+ {X cosiQt + <p x ), Y sm.(Qt + <p v )} (12) 

p{t) = {—vo sin(wi + 0), vo cos(wt + (j))} + 

+ {-Vo, x sin(fit + (p x ), V 0}V cos(Sli + cp y )} , 

(13) 

where X (Yq) is the amplitude, Vo, x = QX (Vo, y = QY ) 
the velocity amplitude and ip x {if y ) the initial phase of 
the macromotion in x (y) direction; 4> is the initial phase 
of the micromotion. The atom starts at rest at the origin, 
hence the initial conditions are p, p = at t = 0. If the 
condition 

uj > n (14) 
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TABLE I: Comparison of numerical results (num) with the 
analytical model (AM); +ab - including refinements (a) and 
(b); +abc - all refinements included 



FIG. 2: (color online) Numerically calculated trajectories in 
the xy plane with the x- and y-positions shown against time 
(left) and parametric plots of the same trajectory in the xy 
plane (middle and right) of a particle initially at rest at the 
origin, after instant switch-on (black lines). The dotted red 
curves correspond to a switch-on time of 3 ps of the TOP field, 
a settling time for the value of Bq as well as the presence of 
gravity. The trap frequencies are uj/2h = 4 kHz and Q/2n = 
394 Hz. Units are scaled to the TOP radius p m . 



is satisfied, the acceleration due to the micromotion dom- 
inates over that of the macromotion. The total accelera- 
tion may be approximated by p ~ F c /m. In other words, 
p points in the direction p m , which is opposite to the di- 
rection of p (as per Fig.[T]). Hence, the initial phase of 
the micromotion is <j) ~ <j) m + tt, where 4> m is fixed by the 
phase of the rotating B m field [26| . Without loss of gener- 
ality we can set 4> m = 0, which means that p m is oriented 
along the positive x direction at t = 0. With this choice 
and setting <p = <p m + w, we find from the initial condi- 
tions: ifxify = 0, Xq = po, and Yq = (l)/Q,)pq. Substi- 
tuting these values in Eq. (I12[) we obtain an equation for 
the macromotion representing an elliptical orbit with its 
major axis oriented perpendicular to the instantaneous 
direction p m of the B m field at t = 0. Since the ampli- 
tude of the macromotion along its major axis is larger 
than the micromotion by the factor u>/Cl, a substantial 
sloshing motion results from the sudden switch-on. Note 
that with increasing oj, the micromotion amplitude po 
decreases like 1/w 2 whereas the amplitude of the slosh- 
ing motion Yq decreases only like 1/oj. For this reason 
the sloshing cannot be neglected in most practical cases 
involving audio-frequency modulation. 



C. Numerical analysis 

To validate the analytical model introduced in Sec- 
tion lll Bl we numerically integrate the classical equations 
of motion in the full time-dependent potential given by 
Eq. |T]) for z = 0, v z = and cp m = 0. The result for 
the trajectory is given in Fig. [2] and exhibits the sloshing 
macromotion described above. The choice of parameters 
is such that it matches the experimental conditions that 
will be presented in Section Hill 

The drawn black lines in Fig. [5] correspond to sudden 
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switch-on of the TOP trap at t = for an atom ini- 
tially at rest at the origin in the absence of gravity. The 
figure clearly shows the micromotion superimposed onto 
the macromotion orientated along the y direction. The 
amplitudes and phases of the macromotion obtained by 
fitting Eq. (|12l) to the results of the numerical calculation 
agree accurately with the analytical model of Section lll Bl 
(see Table IJ) . A more detailed comparison reveals that 
anharmonicities play a minor role; the harmonics of both 
the micro- and macromotion have amplitudes which are 
at least two orders of magnitude smaller than those of 
the fundamentals. 

In order to allow a better comparison with the experi- 
ments to be discussed below we have also performed the 
numerical analysis including several refinements that per- 
tain to our specific experimental situation. These effects 
are: (a) a difference (5y) in gravitational sag between the 
IP and the TOP trap; (b) an exponential switching tran- 
sient of the current in the TOP coils and correspondingly 
in the B m field m/ e = 3 /is); (c) a switching transient 
of - 0.5 ms in the offset field from B = 9.5 x 1(T 5 T at 
the t = to the final value B = 3.1 x 1(T 5 T. 

The initial gravitational sag in the IP trap is 1.2 /im. 
When switching on the TOP, the sag Ay jumps in ~ 3 pis 
to 1.7 /im and settles in ~ 0.5 ms to its final value 1.6 pm. 
due to the decrease of Bq. Thus the gravitational sag in- 
creases jump wise and settles at Sy = 0.4 pm. During 
the same transient the radius of the stationary micromo- 
tion grows from po = 0.21 /tm to po = 0.33 pm and £1 
increases by about 5%. 

The dotted red traces in Fig. [5] correspond to the nu- 
merical calculation including all the above refinements 
relevant to the experiments. We have also investigated 
the effects of gravity, B m -switching and f?o-switching 
separately. We find that the main effect of the settling 
time of Bq is to reduce the amplitude along the major 
axis by ~ 35%. The combined effect of changing gravita- 
tional sag and B m transient is to slightly increase the x 
amplitude as well as to produce a slight tilt angle of the 
trajectory (see right-most panel of Fig. [5]). 

The tilt angle 9 of the macromotion also follows from a 
fit of Eq. (|12p to the numerical results: for known values 
of Xq, Yq, ip x and ip y the angle of rotation ■& to align 
the coordinate system along the major and minor axis is 
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given by 

= | tan" 1 [2sin(^ - ^)X r /(F 2 - A 2 )] (15) 

For ln = the tilt angle equals the rotation angle [9 = 
0). 

The results of a fit of Eq. (fl2|) to the numerical results 
including only the refinements (a) and (b), as well as a 
fit including all three refinements (a), (b) and (c) are also 
given in Table HI Extending the analytical model to in- 
clude the refinements (a) and (b) is straightforward and 
given in detail in Appendix [A] The expressions for the 
amplitudes and phases depend on the model parameter 
To and are given by Eqs. (|A6|) - (|A9|) of the appendix. The 
model parameter r is chosen by ensuring that the value 
of the tilt angle 9 of the model reproduces that of a fit to 
the numerical solution for zero settling time, 9 = 0.0247T. 
This results in To = 3.5 /is. Excellent agreement is ob- 
tained with the numerical model as is shown in Table Q] 
Insight in the cause of the reduction of the major-axis 
amplitude associated with the settling behavior of Bo 
can also be gained using the analytical model. As dis- 
cussed in Appendix [A] the major refinement is change the 
launch speed corresponding to the initially smaller value 
of pq. Although this refinement captures the origin of the 
35% reduction of the major axis amplitude, Table|I]shows 
that the overall agreement with the numerical model is 
less favorable. 



D. Phase jumps 

Let us now analyze how the macromotion can be 
quenched. For a one-dimensional, rapidly-oscillating po- 
tential it was demonstrated in Ref. [6| that the amplitude 
of the macromotion can be quenched by an appropriate 
phase-jump of the modulation field. For the 2D motion 
in a TOP, the success of such an approach is not a priori 
obvious because the phase jumps for the x- and y motion 
cannot be selected independently. Yet, as will be shown 
below, also for the TOP it is possible to quench both 
the Xq- and Yq amplitudes more or less completely by 
imposing a single phase jump A<fi m to the B m field. 

For clarity we first restrict ourselves to the case <j> m = 
and neglect the effects of gravity and switching tran- 
sients. This means that the cloud is launched at t = in 
the vertical y direction with a speed that is equal to Vq, 
the micromotion speed. As can be seen from the trajec- 
tory depicted at the left of Fig. 3 the macromotion speed 
will again be equal to vq when the cloud returns close to 
the origin after an integer number of macromotion half- 
periods. The total velocity p (t) is the vector sum of 
the micro- and macromotion velocities and this quantity 
varies rapidly on a time scale of the micro-motion period. 

The essence of the quenching procedure is to apply 
the phase jump at a time t a chosen in the interval 
t n — At < t <t n + At around times t n — n (7r/fi) corre- 
sponding to a multiple of the macromotion half-period. 
We choose t a such that p (t a ) has a magnitude equal to 




FIG. 3: (color online) Explanatory diagram for the phase 
jump. Left: cloud trajectory (black solid line) along with 
macromotion trajectory (blue dotted line) The black dashed 
lines are the symmetry axes of the trap and the blue arrows 
show the macromotion velocity on crossing the :r-axis. Mid- 
dle: Expanded view of boxed region of the left panel; p (t) 
is the position of the cloud at the time of the phase jump. 
The red dashed (black dot-dashed) circle is micromotion just 
before (after) the phase jump at t = t a . Right: micromotion 
(v) and macromotion(V) velocity vectors add up to the total 
velocity vector p (t) . 



Vq. When the cloud returns at the x axis the micro- and 
macromotion speeds are both vq and hence the resultant 
total velocity can only be equal to vq if the angle between 
the macro- and micromotion directions is either 27r/3 or 
— 27r/3 corresponding to two distinct micromotion phases 
4> a = <f>(t a ) = u)t 2n -i + 4> = ±tt/3 (see Fig. 3-right). In 
other words the micro- and macromotion velocity vec- 
tors form an equilateral triangle. For each of these cases 
a corresponding phase jump exists, A<j> m = ±7r/3 re- 
spectively, such that p m is set perpendicular to p(t a ), 
which sets the macromotion velocity to zero. The result 
is pure micromotion if the orbit into which the particle 
is kicked is centered around the origin. For each of the 
two choices of <p a , pure micromotion results only if the 
macromotion position at the time of the phase jump is 
equal to (ipojO), where the + (— ) sign applies for even 
(odd) n. Complete quenching can be achieved only for 
specific choices of the ratio u>/Cl. The change of orbit 
upon a phase jump is explained pictorially in the middle 
of Fig. 3. 

We now generalize to the case where the ratio w/Cl is 
not precisely fine tuned and allow for the possibility that 
the macromotion speed deviates slightly from the value 
vq assumed above. One can show that, also in this case, 
the maximal reduction in marcromotion energy resulting 
from a phase jump is achieved when the jump is applied 
at a time t a when p(t a ) has a magnitude equal to vq. 
The value of Acf> m is again selected such as to set p m 
perpendicular to p(t a ). By a reasoning similar to the 
case described above we find that the condition of an 
equilateral triangle of the three velocity vectors is now 
replaced by one that is isosceles-triangle condition with 
the micro-motion velocity and p (t a ) both having a mag- 
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nitude vo- This in turn means that the magnitude of the 
phase jump will deviate slightly from the values ±7r/3 
found above. Also, the nearest distance to the x axis 
at which the isosceles-triangle condition can be met is in 
general not equal to zero. This means that some residual 
macromotion will be present after the phase jump, with 
an amplitude given by the distance to the origin of the 
center of the circular orbit into which the cloud is trans- 
ferred by the phase jump. One can show that there is 
always a choice possible where the isosceles-triangle con- 
dition is satisfied such that this distance is approximately 
2p or less. As a consequence, even in the worst case, the 
macromotion amplitude is reduced from (u/Q)po to an 
amplitude of order po- 

The criterion that the acceleration be set perpendicu- 
lar to the total velocity at the time that the macromotion 
speed is equal to w can be expressed by the following 
equation: 



A0„ 



arctan 



Py (*a) 



Px{t a ) 



-l) fe o (16) 



where p x (t a ) = —v sin 4> (t a ) ~ Vb,x sin(Qi Q + ip x ) and 
Py (t a ) = vq cos <j> (t a ) + Vo tV cos(flt a + <p y ) are x- and y 
components of p at time t a and k = 1 for p x (t a ) > 
and k = for p x (t a ) < 0. We return to selection of the 
jump time and the use of Eq. (|16j) when discussing the 
measurement procedure in Section [ill Bl 

Examples of the numerical calculations of the quench- 
ing procedure are shown in Fig. 01 The near complete 
quenching of the macromotion shown in panel (a) is ob- 
tained for 5y = and t = with phase jump A(f> m — 
— 7r/3 at time t a — 3.834 ms in the time interval around 
tj, = In Fig. 0b the refinements (a), (b) and (c) are 

included in the simulation of the experiment. In this case 
the phase jump had to be adjusted to A0 m = — 0.22-7T for 
maximum quenching. Note that the quenching is less 
complete. By adjusting, at constant fi, the micromo- 
tion frequency to u) = 4.068 kHz and the jump time to 
t a = 3.769 ms, complete quenching similar to that shown 
in panel (a) was obtained also when including all refine- 
ments in the numerical model. 



III. EXPERIMENTAL 
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FIG. 4: (color online) Numerically calculated radial trajec- 
tories in the x- (black) and y (red) direction for the same 
trap parameters as used for Fig. [2] with a quenching phase 
jump Acfim applied at optimized t = t a — £3 (three macromo- 
tion half-periods), (a) instant switching, no gravity: A(f> m = 
— 7r/3, t a = 3.834 ms; (b) including switching transients and 
gravity: A(j> m = ~0.22n, t a = 3.834 ms. 



TOP field is produced by two pairs of coils, one in the x 
direction, the other in the y direction as described pre- 
viously in jUJ. The coils consist of only two windings 
to keep the inductance low. The current for the TOP is 
generated by a TTI 4 channel arbitrary waveform genera- 
tor (TGH 1244), amplified by a standard audio-amplifier 
(Yamaha AX-496). The current used is I m = 3.0 A and 
the field produced is B m = 6.8(2) x 10~ 5 T. All mea- 
surements in the TOP are done with 0/27T = 394(4) Hz 
(Bo = 3.1 x 10~ 5 T). Detection is done by time-of-flight 
absorption imaging along the z axis using a one-to-one 
transfer telescope to image the xy plane onto a Princeton 
TE/CCD-512EFT CCD camera with 15 pm pixel reso- 
lution. All measurements are carried out with the same 
flight time AixoF — 23 ms, giving rise to an expanded 
cloud radius of ~ 140 ^m. 



A. Apparatus 

The experiments are done with the apparatus de- 
scribed in detail in [27} and |28| . We produce a BEC 
of 2.5 x 10 5 atoms of 87 Rb in the \F = 2, m F = 2) state 
in a Ioffe-Pritchard trap using radio-frequency (rf) evap- 
orative cooling. The symmetry axis (z axis) of the trap 
lies horizontal with trap frequencies (f2 p /27r = 455(5) Hz, 
fi z /27r = 21 Hz) and the magnetic field offset Bo = 
9.5(3) x 10~ 5 T, a = 3.53 T/m and = 266 T/m 2 . 
The Thomas- Fermi radius of the BEC is 2.2 pm. The 



B. Measurement procedure 

Our experiments on phase-jump-controlled motion in a 
TOP trap are done with the B m field operated at ui/2ir = 
4 kHz. This frequency is sufficiently high (cj/O > 10) to 
satisfy the 'TOP condition' Eg. (fI3]l. The frequency is 
chosen lower than in a typical TOP to ensure that the 
speed of the stationary micromotion, 9 mm/s as esti- 
mated with Eq. (fTT|) . is accurately measurable. In the 
experiments we start with an equilibrium BEC in the IP 
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FIG. 5: (color online) The centroid position after 23 ms TOF 
plotted in camera pixel units against holding time in the TOP 
trap: upper datatset: R x ; lower dataset: R y . The solid lines 
represent the fit of Eqs. (|17[l and (|18Jl to the data. Note that 
by a stroboscopic measurement at 0.25 ms intervals the mi- 
cromotion is eliminated. Each point represents a single mea- 
surement. 



trap described above. At t = we switch on the B m 
field, using Bq to tune the measured trap frequency to 
J7/27T = 394 Hz. As the trap minimum shifts down by 
5y = 0.40 /mi, the initial position of the cloud is slightly 
above the trap center. The 1 /e-switching time of the B m 
field was measured to be t « 3 [is, which corresponds 
to lot « 0.08. When changed, the Bo field settles to a 
new value after a damped oscillation with a frequency 
of 650 Hz and a damping time t' of 0.56 ms. This cor- 
responds to fir' « 0.2. The velocity p of the BEC in 
the radial plane at the time of release is determined by 
time-of-flight absorption imaging along the z axis. For 
the chosen flight time of 23 ms, a speed of I mm/s cor- 
responds to a displacement of 23 /mi with respect to a 
cloud released from the same position at zero velocity. A 
cloud released at rest at time t re i is imaged at position 
Ro = p(t re i) + |p g Ai TOF , where p g is the gravitational 
acceleration. For a finite release velocity p(t re i) the cloud 
will be imaged at R = p(t re i)AiTOF + Ro- 

In practice we may neglect the small variation in the 
release position due to the macromotion, approximating 
p(t re i) ~ p(0), because this variation is smaller than the 
shot-to-shot reproducibility of the cloud position. From 
the model analysis of Section III Bl the variation in re- 
lease position due to the macromotion is estimated to be 
Sp(t re i) < (cu/Q)po « 4 /im. The centroid of the im- 
age of the expanded cloud is determined using a simple 
Gaussian fitting procedure and has a shot-to-shot repro- 
ducibility of ~ 8 /zm, small as compared to the 140 /im 
radius of the expanded cloud. No improvement in shot- 
to-shot reproducibility was found by changing to a higher 
magnification. Since our measurements depend only on 
the position of the cloud center they are insensitive to 
fluctuations in atom number or density. To reconstruct 
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FIG. 6: (color online) Illustration of how to choose the opti- 
mal phase jump and its timing. In both panels: solid black 
curve - experimental conditions; dashed red curve - analyti- 
cal model of Section Til Dl for the case of instant switching - no 
gravity, (a) The total speed of the cloud in units of the mi- 
cromotion speed do (optimal phase jump time t a corresponds 
to p(t a ) — wo); the dashed blue line (scale on right) shows 
the Y component of the macromotion position crossing zero 
at t = t$ (stationary micromotion can be achieved by adjust- 
ing oj such that tz = t a )\ (b) The optimal phase jump as a 
function of jump time as calculated by Eq. (11611 . 



the motion of the condensate in the trap we image the 
cloud at t = ti, where tt is the holding time in the TOP. 
We obtain the release velocity by measuring the x and 
y components of the cloud centroid (R x ,R y ). A typi- 
cal set of data is shown in Fig.[5j The micromotion is 
recognized as the rapid modulation on the slow macro- 
motion. As the frequency of the micromotion is accu- 
rately known we avoid aliasing by sampling the motion 
in steps of 0.025 ms, much shorter than the micromotion 
period. If we wish to look only at the macromotion in a 
stroboscopic manner, we can sample precisely at the mi- 
cromotion period of 0.25 ms, with best results obtained 
when sampling on the crests of the micromotion. Fitting 
the expressions 

R x = — UnAtTOF sin(wi + cf>) — 

-V r o 1 *At T OFBin(nf + v«)+.Ro ia : (17) 
R y = uoAtxoF cos(ujt + 4>) + 

+ Vo <v At T oFCos(m + (py) + Ro,y (18) 

to the data, using the TOP frequency oj and AixoF as 
known parameters, we obtain the amplitudes un, Vo. x , 
Vo t y as well as the macromotion frequency and the 
phases <f>,ip x ,cp y . Note that the fit also yields the refer- 
ence position Ro = {Ro, x , Ra.y} but this information is 
superfluous for the reconstruction of the in-trap motion. 
Once these quantities are determined the motion of the 
condensate in the TOP trap is readily reconstructed with 
Eq. (H2D. 

To investigate the effect of phase jumps, we implement 
the approach described in Section III Dl First we deter- 
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mine for given u> and AixoF all parameters to reconstruct 
the motion with the method just described. This enables 
us to determine the time intervals t n — At < t < t n + At, 
where the cloud returns close to the origin, and choose 
within this interval the time t a , where the total velocity 
p{t a ) has magnitude vq as shown in Fig. [6k,- The red 
dashed lines correspond to the analytical model of Sec- 
tion III Dl for the case of instant switching - no gravity 
(the case of Fig.[4ji). The black solid lines correspond to 
the calculation including all relevant experimental con- 
straints. The phase jump A<f> m that sets p m perpendic- 
ular to p(t a ) is given by Eq. (fTo| . This optimal phase 
jump A(f> m is plotted versus t a in a time interval around 
£3 = 3tt/H in Fig. [6}}. For the case of instant switch- 
ing - no gravity the optimum phase jump is seen to 
be A<p rn = — 7r/3. At the chosen time t a we vary the 
phase jump Acj) m about the value suggested by Eq. (fTB) in 
search for optimal quenching. To reconstruct the residual 
macromotion, we hold the cloud for a variable additional 
time tb, before TOF imaging at time t = t a + %. 



IV. RESULTS AND DISCUSSION 

In this section we show the results obtained with the 
experimental procedure described in the previous section. 
We measured the macromotion induced by switching on 
the B m field for three values of the initial TOP phase, 
<f> m = 0,7r/4, 7r/2. For <f> m — part of the raw data 
are shown in Fig. [5] In Fig. [7] we show the measured 
velocity of the macromotion obtained with the strobo- 
scopic method. The upper and lower panels correspond 
to m = and 7r/4 respectively. The data for cj) m = n/2 
are not shown but are similar to those for <j> — but with 
the roles of x and y interchanged. 

The solid lines in the left panels of Fig. [7] are obtained 
by fitting Eqs. fT7]) and flS]) to the full data including 
micromotion and provide the input for calculating the 
amplitudes. Using the known TOP frequency cj/2tt = 
4 kHz and flight time AixoF = 23 ms, the fit yields 
for the velocity amplitudes, phases, and frequency: vq = 
7.6(2) mm/s, V , x = 0.7(2) mm/s, V , y = 5.6(2) mm/s, 
<j) = 1.00(l)7r, ipx = 0.5(2)tt, tpy = 0.05(2)7T, and n/2n = 
394(4) Hz. The corresponding in-trap amplitudes are 
p = v /u> = 0.30(1) fim, X = V , x /fi = 0.28(7) pm, 
and Y = V , y /fl = 2.3(1) /xm. 

The right panels in Fig. [7] are parametric plots of the 
trajectories obtained by reconstructing the motion in the 
trap from the velocity fits described above. The trajec- 
tories provide a useful way to see the effect of the initial 
phase of the applied B m field and in addition the upper 
panel can be directly compared with the theoretical pre- 
diction shown in Fig. [5] As expected, the orientation of 
the major axis of the macromotion is dependent on the 
initial phase <j) m of the B m field. The small tilt 9 away 
from the direction perpendicular to p m is clearly visible 
and consistent with the calculations for a finite switch- 
on time and the presence of gravity. The value obtained 
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FIG. 7: (color online) The panels on the left show the macro- 
motion velocities (taken with the stroboscopic method) of the 
cloud centroid x- (black) and y- (red) versus time for <f) m = 
and 7r/4. The solid curves are fits of Eqs. (fT7|) and (TTS)) to 
the data. The panels on the right represent the reconstructed 
trajectories in parametric form (in units of the TOP radius 
p m ~ 19.5 [im. The difference in aspect ratio is caused by the 
gravity shift. 



for po is slightly smaller than the value calculated with 
Eq. (|10[) but in view of experimental uncertainties cer- 
tainly consistent with the value of a. The results for 
<p x /ir, (fly/tr, Xq/ pq,Yq/ po and the tilt angle 9 obtained 
for 4> rn = and 7r/4 are given in Table [TT1 For compari- 
son, also the numerical results are included. 

We now turn to the results of a quenching experiment. 
The time t a ~ 3.83 ms and magnitude A<f> m = — 0.227T of 
the phase jump have been chosen to meet the conditions 
necessary to quench the macromotion as introduced in 
Section fill Bl and illustrated in Fig. [6] In Fig. [8] we show 
velocity data taken with the stroboscopic method. For 
t < 3.83 ms the data coincide with those shown in the 
upper panel of Fig. [7] but the solid lines are not a fit but 
represent the macromotion velocity predicted by the nu- 
merical calculation on the basis of the experimental pa- 
rameters. These velocity curves correspond to the macro- 
motion part of the position plot Fig. 2}d and have no ad- 



TABLE II: Experimental results (exp) for macromotion in- 
duced by the switch-on of the TOP field for <j>m = 0,7r/4. 
The data are compared with the results of the numerical cal- 
culation of Section [II CI (num). In all cases the tilt angle has 
be calculated with the aid of Eg . 1151 





4>m 


e/Ti 




Ify/TV 


Xo 1 po 


Y Q /p a 


exp 





0.04(2) 


0.5(2) 


0.05(2) 


0.9(2) 


7.7(3) 


num 





0.016 


0.20 


0.06 


0.82 


6.5 


exp 


tt/4 


0.04(2) 


0.49(2) 


0.12(2) 


6.6(4) 


5.2(3) 


num 


tt/4 


0.013 


0.47 


0.04 


4.95 


4.55 
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FIG. 8: (color online) Measured and calculated velocity of the 
macromotion before and after a phase jump of A<j> m = — 0.227T 
at t a = 3.83 ms for initial phase <f> m = . The open black 
squares (solid red circles) correspond to the measured V x (V y ) 
velocity component (for t a < 3.83 ms the data coincide with 
those of Fig. [7]). Each point represents a single measurement. 
The solid lines correspond to the numerical model without 
any adjustable parameter as described in the text. 



justable parameters. Both experiment and theory show 
pronounced reduction in the amplitude of the macromo- 
tion. Although, the phases of the quenched motion can- 
not be determined convincingly with our signal to noise 
ratio, the agreement between theory and experiment is 
satisfactory. 

In general a jump in the micromotion phase produces 
an abrupt change in macromotion phase and amplitude. 
For the case illustrated in Fig. [5] we obtain a reduction of 
more than a factor of 5 in the amplitude of oscillation in 
the y direction at the expense of only a slight increase of 
the amplitude in the x direction. As a result the macro- 
motion is reduced to the size of the micromotion. The 
energy associated with the macromotion is consequently 
reduced by a factor of about 15, reducing it to a small 
fraction of the micromotion energy. This demonstrates 
that the initial sloshing motion of the cloud can be effi- 
ciently quenched by applying an appropriate phase jump 
angle. As pointed out in the last paragraph of Section 
III Dl we expect that it should be possible to suppress the 
macromotion almost completely by adjusting the micro- 
motion frequency such that t-j =t a . 

Even a small variation in the phase jump magnitude or 
its timing can result in a substantial difference in quench- 
ing efficiency. This is illustrated in Fig. \§\ where we plot 
the ratio of macro- and micromotion energy, 

TP I/ 2 -I- I/ 2 

^macro v 0,x ' v 0,y (-\c\\ 

~E~ ~ v 1 ' [ ' 

as the phase jump is varied in steps of 10 degrees, for 
t a = 1.32 ms and 1.33 ms, where the position and velocity 
criteria are well satisfied. For most phase jumps A</> m 
the result is an increase in energy. The drawn lines are 
the predictions from the numerical model for the same 
conditions at t a = 1.31 ms and t a = 1.32 ms. The plot 
for t a — 1.32 ms shows a deeper reduction than that for 
t a = 1.31 ms, as well as a shifted optimal A</> m . The 
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FIG. 9: (color online) Ratio of macromotion energy over mi- 
cromotion energy following a phase jump plotted against A<f> m 
at 1.32 ms (open black squares) and 1.33 ms (red circles). 
Each data point is obtained from fits as described in the text. 
The horizontal blue dashed line shows the initial value of the 
energy before the phase jump. The solid black lines and dot- 
ted red are the numerical calculation for 1.31 ms and 1.32 ms 
respectively. The inset shows the dependence on jump time 
for fixed value of A0 m . 

common shift of ~ 0.01 ms between the data and the 
numerical results remains unexplained. 



V. SUMMARY AND CONCLUSION 

We have shown that a cold atomic cloud initially at 
rest at the minimum of the effective potential of a TOP 
trap, acquires a macroscopic sloshing motion, in addition 
to near circular micromotion, when the TOP is suddenly 
switched on. The energy associated with this macromo- 
tion is of the same order as the energy of the micromotion 
and the amplitude of the former is larger than that of the 
latter by a factor ~ cj/Cl, We have theoretically described 
the phenomenon and the predictions compare well with 
our experimental results. 

As the micromotion is shared in common mode by all 
trapped atoms, the associated energy does not affect the 
thermodynamics of the cloud in any way. In contrast, 
the macromotion energy is generally unwanted and po- 
tentially harmful. Fortunately, as we have shown, it is 
possible to quench this macromotion almost completely 
and instantly, by applying a suitable and properly timed 
phase jump to the rotating magnetic field that defines 
the TOP. We have shown theoretically that this proce- 
dure works, even for the 2D case of the TOP, which is an 
extension of previous theory describing similar phenom- 
ena in ID |6|, |7| . We have presented a framework which 
allows a deterministic procedure for choosing the optimal 
parameters for the phase jump. Our experiments corrob- 
orate the theoretical model for the TOP in a quantitative 
manner. 

The macromotion induced by the switch-on and the 
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subsequent possibility to alter this motion by phase 
jumps have several consequences, some of which we now 
briefly mention. For example, the sloshing motion may 
affect the time of flight imaging once the fields have been 
switched off. When comparing TOF-images for different 
holding times it is in general not sufficient to synchro- 
nize the release time to the micromotion period. The po- 
sition after TOF can be easily polluted by the non-zero 
macromotion, which evolves asynchronously with the mi- 
cromotion. The time scales in this experiment are on the 
order of a few macromotion periods. The physics of in- 
terest of the cloud is usually seen on much longer time 
scales of hundreds of such periods. On these longer time 
scales, the presence of even small anharmonicities can 
lead to the conversion of macromotion energy into heat. 
The macromotion can be of order of the chemical poten- 
tial which can have consequences for the stability of the 
condensate. 

The possibility to excite or quench macromotion by 
phase jumps of the rotating field is a valuable feature of 
the TOP trap that has received little attention in the 
literature. Our work shows that this feature is well un- 
derstood and can be applied in a well-controlled man- 
ner. We have primarily focussed on quenching with a 
single phase jump. However, the reverse effect in which 
the macromotion is excited may prove equally useful in 
some experiments. Also the consequences for multiple 
phase-jump applications deserve attention in this respect. 
We established numerically that it should be possible to 
excite or deexcite large macromotion with a series of 7r 
phase jumps at intervals of the macromotion half-period. 
At each of these phase jumps, either component of the 
macromotion velocity can be increased or decreased by 
~ 2vq. Being outside the primary focus of this paper, we 
do not further elaborate on this interesting topic. 
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Appendix A: Analytic Model 

For arbitrary 4> m the position p' — {x' , y'} with respect 
to a coordinate system rotated over an angle <p m is 

x'{t) = -p cos tut + X„ cos(fM + ip' x ) (Al) 
y'(t) = -p sin + F ' sm(Qt + ip' y ), (A2) 

where X' Q , Yq are the amplitudes and Lp' x ,ip' y the phases 
with respect to the rotated axes. Taking the time deriva- 



tive and using the initial conditions p',p = at t = 0, 
yields: cp' x = ip' y = 0, X = p , Y ' = (w/O)p . This 
corresponds to an ellipse with its major axis oriented 
perpendicular to the instantaneous direction p m = p m 
of the B m field at t — 0. 

For exponential switch-on of the B m field with 1 /e time 
To, we have for the acceleration in the 'primed' coordinate 
system 

= cj 2 (1 — cxp[— i/r ])po cosut — 

-f2 2 A cos(f2i + <^) (A3) 
y'{t) = cj 2 (1 — cxp[— £/to])po sinwi — 

-n 2 Y^sm{nt + ^ y ), (A4) 

where, during switch-on, X'q, Yq, ip' x , ip' y and Q are func- 
tions of time. Since f! -C w we may approximate, for 

Ujt < 1, 

x'(t) ~ w 2 (l - exp[-i/To])po coscot > y'(t). (A5) 

This shows that the switch-on profile mainly affects the 
acceleration in the x' direction because this is the initial 
direction of acceleration. By the time sinwi is sufficiently 
large to make y' non-negligible, the switch-on transient is 
already finished. For ut C 1, the velocity in the x' direc- 
tion is given by x'(t) ~ uj 2 (t — t + t exp[— t/r])p . This 
expression suggests to approximate the switch-on profile 
by a step function at t = To ~ t and to treat X' Q , Yq, ip' x , 
ip'y as constants for t > tq. This 'delayed sudden-step ap- 
proximation' is equivalent to imposing the boundary con- 
ditions p' ,p = at t = To. In this approximation we ob- 
tain for the phases and amplitudes ip' x — tan _1 [w 2 To/r2], 
V ' y = 0, X'q = p (l + ^rH^fl\ Yq' - (w/n)po for 
t > tq. The phase development lutq due to rotation of 
the B m field during switch-on will appear as a rotation 
of the major and minor axes of the macromotion with re- 
spect to the primed coordinate system defined by t = 0. 
The optimal value for To can be determined by comparing 
the predictions of the analytical model with the results 
of a numerical calculation (see Section ["II C[) . 

In the presence of gravity, the above analysis re- 
mains valid as long as the radial frequency does not 
change substantially during switch-on of the B m field; 
in principle f2 can be kept constant by simultaneously 
switching Bq and B m in such a way that the quantity 
B /(l + ^B^/Bq) equals the value of B before the B m 
field was switched on. In case of a small and fast change 
in the above model can be adapted by changing the 
initial conditions to p — {0,Sy},p — at t — To, where 
Sy = Ai/top — Aj/iP is the difference in gravitational 
sag. Using the adapted boundary conditions we obtain 
in the limits {ga/ Bq) 1 / 2 < < w and ujtq -C 1 for the 
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amplitudes and phases of the macromotion 
X Q = po[l + {uj 2 /Q 2 - 1) sin 2 ( W r + m )] 1/2 (A6) 

*0 = P [l + (w 2 /0 2 - 1) COS 2 (WT + <j)m) 

+ Sy 2 /p 2 + 2(6y/p ) sin( W r + 0™)] 1 / 2 (A7) 

ip x = —£It + tan -1 [(w/f2) tan(wr + <fi m )] + nir (A8) 
ip y = —Qtq + tail [(fi/w){tan(u;To + </> m ) 

+ (Sy/po)/ cos(wt + </> m )}] + rwr. (A9) 

where n — for |wro + m < 7r /2, n = 1 for |wro + <^ m | > 
7r/2. For </> m = 0, these equations coincide with the equa- 
tions for Xq, Yq, ip' x and ip' y in the primed coordinate sys- 
tem. Analyzing the limit ujtq — > for the case Sy ~ po 
(typical for our experimental conditions) we find Xq ~ po 
and Yq ~ (ui/Q)po for </> m = and Ao ~ (w/Q)po and 
y — 2po for m = 7r/2. Thus, we deduce that gravity 
can have a substantial influence on the amplitude of the 
macromotion along its minor axis but not on the am- 
plitude along the major axis. For known values of Xq, 
Yq , (p x and ip y the angle of rotation $ to align the coor- 



dinate system along the major and minor axis is given 
by Eq. (fT5j) . For Ht <C wt <C 1 and in the absence 
of gravity (5y = 0) the delayed sudden-step gives rise to 
a small rotation Aip ~ lutq, independent of (f> m . For 
Sy ~ po gravity gives rise to an additional contribu- 
tion to this rotation, which is minimal for <p m = ir/2, 
where the macromotion is launched perpendicular to the 
gravity direction and Eq. (|15|) can be approximated by 
<p ~ ojtq[\ + (1 + Sy/p ) (fi/w) 2 ]. The contribution is 
maximal for <p m = 0, where Eq. ()15|1 can be approximated 
by Lp ~ wr - (tt/uj) 2 (Sy/po) 2 (1 + ut ). 

Insight in the dependence of A and lo on the settling 
behavior of Bq can be obtained from Eq. (|A5[) , which 
shows that the initial acceleration and, hence, the launch 
speed scales with the initial value of po. Therefore, most 
of the settling behavior is captured by using the initial 
value of po in Eqs. (|A"6|) and (|A7|) . The 5% change in tt 
requires a further refinement of the model. This cannot 
be implemented without sacrificing the simplicity of the 
model and is not pursued here. 
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